Granular superconductivity in the cuprates evinced by finite size effects in the specific 

heat and London penetration depths 
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We review and refine the finite size scaling analysis of specific heat and London penetration 
depths data of cuprate superconductors and compare it to the analysis of specific heat measure- 
ments near the superfluid transition of 4 He confined to l^im 3 cylindrical boxes. This system crosses 
from 3D to OD behavior near the transition. This has a marked effect on the specific heat as seen 
by a pronounced rounding of the maximum and a shift to a temperature lower than the transi- 
tion temperature of the bulk system. The region in between the 3D to OD crossover uncovers the 
contributions from the surface and the edges of the cylindrical boxes. Our finite size scaling anal- 
ysis of the specific heat and London penetration depth of high quality YBa2Cu307_f and Bi2Sr2 
CaCu2C>8+<5 single crystals uncover essentially the same crossover phenomena, including evidence for 
surface and edge contributions. This implies that the bulk samples break into nearly homogeneous 
O ' superconducting grains of rather unique extent and with that granular superconductivity. 
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Since the discovery of superconductivity in cuprates by Bednorz and Muller [1] a tremendous amount of work has 
been devoted to their characterization. The issues of inhomogeneities, granularity, and their characterization are es- 
sential for several reasons, including: (I) If inhomogeneity is an intrinsic property, a re-interpretation of experiments, 
measuring an average of the electronic properties, is unavoidable. (II) Inhomogeneity may point to a microscopic 
phase separation, i.e. superconducting grains, embedded in a non-superconducting matrix. (Ill) There is neutron 
' spectroscopic evidence for nanoscale cluster formation and percolative superconductivity in various cuprates [2,3]. (IV) 
1—1 , Nanoscale spatial variations in the electronic characteristics have been observed on the surface of Bi2Sr2CaCu208+i5 
with scanning tunnelling microscopy (STM) [4-7]. They reveal a spatial segregation of the electronic structure into 
, 3nm diameter superconducting domains in an electronically distinct background. (V) Remarkably, in the experi- 
■^j- ■ mentally accessible temperature regime, the critical properties of the superconductor to normal state transition are 
s * , found to fall into the 3D-XY universality class, except from a rounding of the transition [9,8,10-25]. This, however, 
is the result expected according to the Harris criterion [26] , which states that short-range correlated and uncorrelated 
disorder is irrelevant at the unperturbed critical point, provided that the specific heat exponent a is negative. Since 
in the 3D-XY universality class a is negative [27], the rounding of the transition appears then to uncover a finite 
size effect [28,29]. Indeed, consistency of the rounded transition in the specific heat and penetration depths with a 
■ finite size effect has been established in a variety of cuprate superconductors. The finite size scaling analysis revealed 
that homogeneity is restricted to grains with nanoscale extension [22,30-34]. (VI) On the other hand, single crystals 
of cuprate superconductors conduct in the normal state along any direction. Accordingly, bulk superconductivity 
requires the grains to percolate. (VII) If this finite size scenario holds true, the finite size behavior of the supercon- 
ductor to normal state transition should resemble features of the superfluid transition of confined He, which falls 
into the 3D-XY universality class, as well [35-38]. In particular, the finite size effect should also provide information 
O ■ on the geometry of the grains, e.g. their surface, edges, corners, etc. (VIII) Noting that in Yi- x Pr x C\i30j-s [30] and 
YBa 2 Cu40g [34] the transition temperature increases with reduced extent of the grains along the c-axis, the geometry 
turns out to be a relevant feature. Indeed, if the grain is in contact with a superconducting layer with a higher "T c ", 
superconductivity is predicted to be enhanced and T c increases with reduced radius of curvature of the grains [39] . 

In this paper we review and refine the evidence for finite size behavior in the specific heat and penetration depth 
data of cuprate superconductors. Since 4 He confined to uniform small dimensions has been the subject of fundamental 
research for over three decades, we sketch next the finite size scaling analysis of the specific heat of 4 He confined to 
cylindrical boxes whose diameter and height are 1 fim [38] . This provides information on a system where confinement 
is made by design and modifies its critical behavior, while other effects such as disorder are not present. Close to the 
bulk transition temperature T c the singular part of the specific heat adopts the form 

c b (t) = c bs (t) + B ± , c bs (t) = — \tl a . (1) 
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t = 1 — T/T c is the reduced temperature, A* 1 the critical amplitude and a the critical exponent of the specific heat, 
while B arises from the background, and ± = sgn(t). 

Since superconductors in the experimentally accessible critical regime and 4 He fall into the 3D-XY universality 
class, with known critical exponents and critical amplitude combinations, we take these properties for granted [27]. 
They include the exponents 

a = 2 - iv = -0.013, v = 0.671, (2) 

and the critical amplitude combinations 

A±vt = {R±)\vt = {£)\ (3) 

where 

A + 

— = 1.07, Br = 0.815, R + = 0.361. (4) 

£q is the critical amplitude of the correlation length, diverging in the bulk system as ^ = ^ \t\~ u , and — (C^) 3 
is the correlation volume. 
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FIG. 1. (a) Specific heat c of 4 He confined to 1/im 3 cylindrical boxes versus |t|compared to the bulk data. The lines represent 
the bulk data and the open symbols are the confined system. Taken from Kimball el al. [38]. The insertion shows c versus t 



around t„ 



-3.9 10 . (b) Corresponding finite size scaling function g (y) — 1 versus y = |t/£ p |. The red curve indicates the 



limiting behavior g+ (\y\ -» 0) = g£ |yp in terms of g + = 0.993 \y\ ■ The black line, g + (y) 
the surface and the blue one, g + (y) — 1 = 0.04 |j/|~ 4,/3 the edge contributions. 
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Fig. la shows the specific heat of 4 He confined to 1/oti 3 cylindrical boxes and contrasts it to the unconfined system, 
the solid lines [38]. The lower branch corresponds to T > T c and the upper to T < T c . Notice at large values of \t\, 
the data match the bulk data. As one moves closer to the bulk transition, the confined systems specific heat begins 
to systematically deviate from the bulk data. This is predicted by finite-size scaling theory [29,28]. As the system 
feels its finite size when the correlation length £ becomes of the order of the smallest confining length, the specific 
heat adopts the scaling form, asserting that 



c fc (t)~c(t,L) 
\c bs (t)\ 



g ± (y)-l, y = t/\t p \, \t p \ = fa/L) 



(5) 



<7 ± (y) is the finite size scaling function which satisfies: g (—1) = 1 at T = T p , where £ (t p ) = £ |t p | v = L, so that 



y = y P 



T, while g (y —> 0) 



of the confined system. g + (y). 
indicates the 0D limiting behavior g + (\y\ — > 0) 
smaller than L, so that \t\ 
D = 3as [29] 



\y\ holds for a < and g ± (\y\ 
derived from the data of Kimball et al. 



y\ 



oo) = 1. L is the smallest spatial length 
3], is displayed in Fig. lb. The red curve 
In between these limits, where the correlation length is 



> \t p \, the free energy density (and other thermodynamic quantities) can be expanded in 



/ (t, L) - h (t) = -f s (t) + -5/ a (t) + -,f c (t) 



L 2 ' 



L 3 ' 



(6) 
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The contributions / s , f e and f c can be attributed to the surfaces, edges and corners of the confinement, ft is expected 
that each one of these terms contribute in a limited region of the scaling variable. Since f s and f e scale as 
and their contribution to the specific heat c oc —d 2 f/dt 2 scales as |i|~ Q ~ I/ and \t\~ a 2v , respectively. Hence 
the surface term should behave as \t/t p \~ L ' and the edge term as \t/t p \~~ v in the scaling form, Eq. (5). A glance 
to Fig. lb shows that as one proceeds to smaller values of the scaling variable y = \t/t p \ there is clear evidence 
for crossovers from the edges, oc It/tp] -4 , to the surface, oc \t/t p \~ 2 ^ , contribution and finally to the asymptotic 
behavior, g + (|y| — > 0) = g$ |y|"' a '- Thus, the finite size scaling function does not uncover the crossover from 3D to 
0D behavior only, but allows to identify surface and edge contributions, as well. 

We have seen that 4 He confined to 1/j.m 3 cylindrical boxes crosses over from a 3D to a 0D behavior near the bulk 
transition. This has a marked effect on the specific heat as seen by the pronounced rounding of the peak and its 
shift to a temperature lower than the transition of the bulk system. Furthermore, the region in between the 3-D 
to 0D crossover uncovers the contributions from the surface and the edges of the cylindrical boxes. Since disorder 
is not present, the modifications of the critical behavior are entirely due to the confinement. Analogous behavior 
in cuprate superconductors would imply that the bulk breaks into homogeneous superconducting grains with rather 
unique extent, embedded in a nonsuperconducting matrix. 

In Fig. 2 we show the specific heat data of Charalambous et al. [40] of a high quality YI^OusC^-^ single crystal 
in terms of the specific heat coefficient c/T versus \t\. Comparing Figs. la and 2 we observe: (i) in YBa2Cu307_^ 
the rounding of the transition sets in at much higher reduced temperatures; (ii) otherwise the rounding is remarkably 
analogous to that of 4 He confined to l/im 3 cylindrical boxes. Noting that in 4 He and optimally doped YBaaCuaCv-is 
the critical amplitude ^ are comparable, the relation \t p \ — /L) 1 ^ (Eq.2) indicates that the difference in the 
onset of the rounding is due to a limiting length scale, considerably smaller than the lfim in 4 He. 
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FIG. 2. Specific heat coefficient c/T of the high quality YBa2Cu307_<s single crystal (sample UBC2) versus \t\ taken from 
Charalambous et al. [40] . The solid lines represent the expected behavior of a fictitious homogeneous bulk sample according to 
Eq.(l). The insertion shows c/T versus t around t p ~ —5.8 1CF 4 . 



To substantiate and refine the consistency with a finite size effect further we displayed in Fig. 3 the scaling functions 
g~ (y) — 1 and g + (y) — 1. As one proceeds to smaller values of the scaling variable y = \t/t p \ there is, in analogy 



to 4 He confined to 1/im 3 cylindrical boxes (Fig. lb), clear evidence for crossovers from the edge, oc \t/t p 



-4/3 



to the 



surface, oc \t/t } 



,-2/3 



contributions and finally to the 0D behavior, g + (\y\ — > 0) 
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1 versus y = t/t p derived from the data shown in Fig. 2. The red curve indicates the limiting behavior 
"' in terms of g + = |i/| _0 ' 013 . The black line, g~ (y) — 1 = 0.01 |y| -2 ^ 3 uncovers the surface and the blue 
0.016 jyj~ 4//3 the edge contributions. (b) g~ (y) — 1 versus y — \t/t p \ derived from the data shown in Fig. 2. The 



FIG. 3. (a) g~ (y) 

g~ o) = go \y 

one, g~ (y) - 

red curve indicates the limiting 0D behavior, g + (\y\ — > 0) = |y| -0 ' 013 . The black line, gr" 1 
the surface and the blue one, g + (y) — 1 = 0.0016 \y\~ 4 ^ 3 , to the edge contributions. 
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corresponds to 



When this finite size scenario holds true, the London penetration depth should not only exhibit a rounded transition, 
consistent with a finite size effect [30- 34], but uncover the surface and edge contributions, as well. Considering again 
the 3D-XY critical point, extended to the anisotropic case, the penetration depths and transverse correlation lengths 
in directions i and j are universally related by [22,41] 
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Wir 3 k B T 



A* (T) A, (T) $ 2 



v %'(/)c; (/•)' 



(7) 



< L\, where 



A limiting length scale Li in direction i implies that does not diverge but is bounded by = )' 

i^j^fc. £| denotes the transverse and £^ = the corresponding real space correlation length. The resulting 

finite size effect has been established in a variety of cuprates [30-34]. Here we refine these studies to identify surface 
and edge contributions. As first example we consider the high-quality Bi2Sr2CaCu20s+5 single crystal data of Jacobs 
et al. [42] shown in Fig. 4. The solid curve indicates the leading critical behavior of the fictitious homogeneous system 
with T c — 87. 5K, while the rounded transition points to a finite size effect. A characteristic feature of a finite size 
effect in 1/A^ fc (T) is the occurrence of an inflection point giving rise to an extremum in d {\ 2 ab (T = 0) /\ 2 ab (T)) jdT 
at T p . Here Eq.(7) reduces to 
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16TT 3 k B T p 



(T p ) A a (T p ) A fc (T p ) 



(8) 



The data shown in Fig. 4 exhibits in d {\ 2 ab (T — 0) /X 2 ab (T)) JdT at T p w 87K an extremum and with that an inflection 
point in l/\ 2 ab (T). With \ ab (T = 0) = 1800A as obtained from /iSR measurements [43] and \ 2 ab (T = 0) /X 2 ab (T P ) = 
0.066 we find L r w 68A. 
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FIG. 4. (a)Microwave surface impedance data for X 2 ab (T = 0) /X 2 ab (T) (O) and d (X 2 ab (T = 0) /X 2 ab (T)) /dT (□) versus T of a 
high-quality Bi 2 Sr 2 CaCu2 08+5 single crystal taken from Jacobs et al. [42]. The solid line is X 2 ab (0) /X 2 b (T) = 1.2 (1 - T/T c ) 2/3 
with T c — 87. 5K and the dashed line its derivative indicating the leading critical behavior of the fictitious homogeneous system. 
The rounded transition exhibits an inflection point at T p » 87K, where d (A^ b (0) /A^ b (T)J /dT is maximum, (b) Scaling 
function g (y) = (Xo ab /X ab (T)) 2 \t\~ v versus y = tj \t p \ for the data shown in Fig. 4a. The solid line in the insertion is Eq.(10) 
with go = 1.6. It indicates the flow to the 0D behavior. 



To strengthen the evidence for a finite size effect we explore next the consistency of the data with the corresponding 
finite size scaling function. In the present case it is defined as 

For t small and L c — > oo, so that ±y — > oo it should tend to g~ (y — > — oo) = 1 and g + (y — > oo) = 0, respectively, 
while for t = and L ^ it diverges as 

g ± {y^0) = gt\yr ■ (10) 

In this 0D limit, (A 0a h/A a h (T c , L c )) 2 = g^ \t p \ v = g £,Q C /L c . Moreover at t p , y = -1 and d (X 0a b/Kb (T, L)f fdt = 0. 
The scaling function shown in Fig. 4b, derived from the data shown in Fig. 4a, is apparently fully consistent with the 
limiting behavior of the finite size scaling function and uncovers, in analogy 4 He, confined to cylindrical boxes, a 3D 
to 0D crossover. On the other hand, in analogy to the specific heat, refined information may be attained in between 
these limits. Here the correlation length is smaller than L c , |t| > \t p \, and the contributions from the surface and edges 
should become observable. Since f s and f e scale as 1/ (l<£ 2 ) and 1/ (£ 2 £), their contributions to 1/A 2 oc / £ 2 scale 
as l/L and £/L, respectively. Hence the surface term should behave as \t/t p \~ v and the edge term as \t/t p \~ 211 in the 
scaling form, Eq. (9). A glance to Fig. 5 shows that as one proceeds to smaller values of the scaling variable y — \t/t p \ 
there is indeed evidence for crossovers from the edges, oc \t/t p \~ 4 ^ 3 , to the surface, oc |t/t p p 2 ^ 3 , contributions and 
finally to the 0D behavior, g + (\y\ — > 0) = g$ |yp 2 ^ 3 . 



(9) 
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FIG. 5. (a) Scaling function g + (y) versus y = \t/t p \ derived from the data shown in Fig. 4. The blue line, g + (y) = 1.7y -4 ' 3 
corresponds to the edge and the black one, g + (y) = 1.6y~ 2 / 3 , to the surface contribution, (b) g~ (y) — l versus y — \t/t p \ derived 
from the data shown in Fig. 4. The blue line, g~ (y) — 1 = 0.7y~ 4 ^ 3 , indicates the edge, the black one, g~ (j/) — 1 = 1.2j/ _2//3 , 
the surface contribution and the red curve indicates the flow to the OD behavior, g~ (y) = 1.6j/ _2//3 . 

As second example we treat the in-plane penetration depth data of Kamal et al. [14] for a high-quality YBa 2 Cu 3 06.95 
single crystal. In Fig. 6a we displayed the data in terms of 1/A^ fc (T) and d (l/A„ 6 (T)) jdT versus T. The red and 
blue curves indicate the leading critical behavior of the fictitious homogeneous system with T c = 88.65K, while the 
rounded transition points to a finite size effect. The occurrence of an extremum in d (l/A^ 6 (T)J jdT at T p » 88.52K 
is a characteristic feature of a finite size effect in 1/A;? 6 (T), giving rise to an inflection point in 1/A 2 ^ (T) at T p . Using 
1/A^ b (T p ) — 0.918 (fim)~ 2 , we deduce with the aid of Eq.(8) for the limiting length scale along the c-axis the estimate 
L c = 154A. The associated scaling function g(y), defined in Eq.(9), is shown in Fig. 6b. As one proceeds to smaller 
values of the scaling variable y there is in the lower branch (T > T c ) clear evidence for a crossover from the edge, oc 
\t/t p \ 4 ^ 3 , to the surface, oc \t/t p \ 2 ^ 3 , contributions. 




FIG. 6. (a) In-plane penetration depth data for 1/A 2 b (T) (O) and d (l/A 2 b (T)) /dT (□) versus T of a high-quality 
YBa 2 Cu 3 06.95 single crystal taken from Kamal et al. [14]. The solid line is l/X 2 ab (T) = 68(1 - T/T c ) 2/3 with T c = 88.65K 
and the dashed line its derivative indicating the leading critical behavior of the fictitious homogeneous bulk system. The 
rounded transition exhibits an inflection point at T p ~ 88.52K, where d (l/A 2 6 (T)) /dT is maximum, (b) Scaling function 
g (y) — (Xoab/Xab (T))' 2 \t\~ v versus y = t/ \t p \ for the data shown in Fig. 4a. The black line, g + (y) = 0.36y _2//3 , indicates the 
surface and the blue one, g + (y) = 0.28y _4,/3 , the edge contributions above T c . 



We reviewed and refined the finite size scaling analysis of specific heat and London penetration depths data of cuprate 
superconductors and compared them to the analysis of specific heat measurements near the superfluid transition of 
4 He confined to 1/xm 3 cylindrical boxes. The 4 He system provides information on a system where confinement is made 
deliberately. It modifies its critical behavior while other effects such as disorder are not present. We have seen that 
this system crosses from 3D to OD behavior near the transition. As a result of this crossover, the specific heat peak 
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exhibits a pronounced rounding and the maximum shifts to a temperature lower than the transition temperature 
of the bulk system. The region in between the 3D to OD crossover uncovered contributions from the surface and 
the edges of the cylindrical boxes. The analysis of the specific heat and London penetration depth of high quality 
YBa2Cu307_,5 and I^S^CaC^Os+i single crystals uncovered essentially the same crossover phenomena, including 
evidence for surface and edge contributions. No new critical behavior emerges, but 3D-XY scaling is achieved with 
the variable £/L, where L is the smallest confining dimension. This implies that the bulk samples break into nearly 
homogeneous superconducting grains of rather unique extent and with that granular superconductivity. 
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